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Bekenstein's alternative to general relativity, TeVeS, reduces to Modified Newtonian Dynamics 
(MOND) in the galactic limit. On cosmological scales, the (potential well<->overdensity) relationship 
is quite different than in standard general relativity. Here we investigate the possibility of cross- 
correlating galaxies with the cosmic microwave background (CMB) to probe this relationship. At 
redshifts of order 2, the sign of the CMB-galaxy correlation differs in TeVeS from that in general 
relativity. We show that this effect is detectable and hence can serve as a powerful discriminator of 
these two models of gravity. 

I. INTRODUCTION 

The standard theory of gravity, general relativity (GR), coupled with the standard model of particle physics cannot 
account for many cosmological observations, ranging from the dynamics of stars in galaxies to the large scale expansion 
of the universe. This discrepency may well be resolved by adding to the standard model of particle physics. For 
example, dark matter in the form of supersymmetric particles could account for flat rotation curves and the large 
scale structure in the universe. Dark energy in the form of a cosmological constant or a scalar field could explain the 
acceleration of the universe. This scenario is commonly called the "ACDM model" . 

Alternatively, general relativity may be wrong. Perhaps the theory needs to be modified on large scales to account 
for cosmological observations. Many attempts have been made along these lines [1-7], to explain either (or sometimes 
both) the acceleration of the universe and the phenomena usually associated with dark matter. One such attempt is 
MOdified Newtonian Dynamics (MOND) initially proposed [8] to explain flat rotation curves in galaxies. Bekenstein [9] 
has recently introduced TeVeS, a theory which reduces to MOND in the appropriate limit but is robust enough to 
make predictions in many different arenas. 

Several groups [10-13] have now studied the evolution of cosmic structure in TeVeS and the theory seems to be 
holding up. In some senses, this is surprising [14] since any no-dark matter theory faces the significant hurdle of 
explaining why the density field has gone nonlinear. Observations of the cosmic microwave background (CMB) pin 
down the baryon inhomogeneities at recombination to one part in 10 4 — 10 5 , and in standard general relativity, these 
inhomogeneities grow by only a factor of a thousand from that epoch until today. 

As shown in [12], overdensities in TeVeS grow faster than in standard general relativity (GR). The mechanics of 
this is subtle: the vector field in TeVeS develops an instability and this growth in turn sources a difference between 
the two Newtonian scalar potentials. The two potentials $ and ^ characterize scalar perturbations to the metric: 

ds 2 = -dt 2 [1 + 2*] + a 2 (t)dx 2 [1 + 2$] . (1) 

In ACDM with standard GR, these two potentials are nearly equal to one another, so the TeVeS situation of unequal 
potentials is unfamiliar territory. In particular, the relation between gravitational potential and overdensity is quite 
different than in standard GR. Fig. 1 illustrates this. 

In this paper, we demonstrate that one can test this prediction of TeVeS (for a different test, see [15]) by cross- 
correlating a deep galaxy survey with the CMB. A number of groups have already succeeded in cross-correlating tracers 
of the density field with the CMB [16-19]. The non-zero correlation is a result of the time-dependent gravitational 
potentials at late times because the universe is not completely matter dominated. In TeVeS, the potentials are 
determined not only by the matter overdensities but also by the behavior of the vector field. So the pattern of the 
cross-correlation changes. Most strikingly, at redshift of order 2—3 the cross-correlation changes sign. Observations 
of such a sign change would be a smoking gun signature of a modification of gravity. 

Before proceeding with the calculation, we note that we focus on TeVeS as a fairly well-defined test model. 
Bertschinger [20] has pointed out however that generic modifications of gravity lead to changes of the sort discussed 
here. So we expect the cross-correlations envisioned here to be powerful tools cosmologists will use to discriminate 
(general relativity + dark matter + dark energy) from many alternative theories of gravity. This idea has already 
been applied in the context of DGP and f(R) gravity for which specific CMB-galaxy cross-correlation signatures at 
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FIG. 1: The "poisson ratio" — fc 2 ($ + *5f) / (Cl^Si, + Q V 5 V ) for the two TeVeS models considered in section III. In these units, the 
Poisson ratio is constant and equal to 3 in standard ACDM cosmology (solid line). This is true in TeVeS only for the largest 
scales, where the standard term in the Poisson equation dominates. On smaller scales the vector field perturbations dominate 
and the Poisson ratio exhibits a much more complicated behavior (see §111 for more details). 



high rcdshifts have been predicted by the authors of [21, 22]. Recently, more general parametrizations of modifications 
to gravity have been proposed [23, 24]. These authors concentrated on modifications of gravity intended to explain 
late-time acceleration. However, TeVeS, with the aim of explaining the breadth of cosmological observations without 
dark matter, shows much more severe deviations from general relativity. In particular, the post-Newtonian parameter 
introduced in [24] is applicable if it is close to unity, which is not the case for TeVeS throughout the whole history of 
the universe. Hence, it is preferable in this case to compare the theory directly with observations. 



II. INTEGRATED SACHS- WOLFE EFFECT 



A CMB photon traveling through time- varying gravitational potentials $ and \& in the metric (1) will suffer a change 
in energy, resulting in a change in the observed CMB temperature: the so-called Integrated Sachs- Wolfe (ISW) effect 
([25-27], see [28] for a review). The fractional change in CMB temperature in the direction n is given by an integral 
along the line of sight back to the redshift of last scattering : 

~T = I dZ d~z (<1>(r ' Z) + * (r ' Z)) |r=x(z)fl (2) 

where x( z ) is the comoving distance out to redshift z. The evolution of the gravitational potentials therefore determines 
the ISW effect. On the other hand, the fractional overdensity of galaxies along the line of sight is 

/>oo 

S g (&)= / dzW(z)6 g ( X (z)n;z) (3) 
Jo 

where W(z) is the redshift distribution of galaxies (normalized so it integrates to unity), determined both by the 
intrinsic redshift distribution of galaxies and by the properties of the survey at hand. 

Since we expect an overdensity of galaxies to be correlated with potential wells, we expect a non-zero correlation 
between the observed temperature anisotropy and the observed angular galaxy distribution. This correlation is 
however model-dependent and can be used to distinguish between different theories of gravity. 
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If we focus on large scales, the fractional galaxy overdensity in Fourier space is 

5 g (k,z) = b(z)S(k,z = 0)D(k,z) (4) 

where b(z) is the time-dependent but scale-independent bias, and D(k,z) is the linear growth function, normalized 
to unity today. In standard ACDM, the growth function depends on z only, but more generally (and in TeVeS in 
particular) it can have a scale dependence. With this restriction to large angular scales, the angular cross-power 
spectrum at a multipolc I is 

2 r°° 

C gT (l) = - dkk 2 P(k)6 9 (l,k)S lsw (l,k) (5) 



7T ^ 

where P(k) is the matter power spectrum today, and the two weighting functions for galaxies and the ISW effect are: 

/>oo 

S g (l,k) = / dzW(z)b(z)D(z,k)ji(k X (z)) (6) 
Jo 

<fisw(i,*0 = J o dz-{- g^Jy ' yMkxiz))^ J o dzD lsw (k,z)Mk X (z)) (7) 

Here, S g (l, k) gives the contribution of the modes with wave number k to the projected overdensity of galaxies for the 
spherical harmonic I and ji is the spherical Bessel function. 

On small angular scales, i.e. large I > 10, the Limber approximation can be used: 



/ 



dk k 2 f(k) ji(k X )ji(k X ') ^^\fik = {l + 1/2) I x) <5 D irac(x - X% (8) 

2 X 



where f(k) is a slowly varying function. This simplifies the expression for the galaxy-CMB cross-power coefficients to 
a single integral (we set c = 1): 

<V(0 - J dz^D(z, k ± ) b(z) W(z) D ISW (k ± ,z) P{k ± ), k ± = l ~^ L - (9) 

Once C 9 t(1) is calculated, the galaxy-CMB angular correlation function w 9 t(0) is given by (e.g., [29]): 

00 o/ 4- 1 

^^) = E^^w (io) 

We use the exact expression (equation (5) and following) up to I = 10. For larger I, we use the Limber approximation, 
as the differences to the exact calculation are very small, about 10~ 3 , as we explicitly verified. 

Equation (5) and following, or cquivalently equation (9), state that the C 9 t(1) are determined by a weighted 
integral of £>isw(&, z) along the line k = (I + l/2)/x(z) in the (k, z)-plane. The function Z?isw defined in equation (7) 
represents the response of the potentials to the growth of matter overdensities. In the case of standard ACDM, it is 
given by the Poisson equation: 

fc 2 $= i -Hln m {l + z)5{k,z), and*-* D isw (k,z) = -^H 2 fl m ^(1 + z)D(z) (11) 

The left panel of Fig. 2 shows Disw as a function of k and z in the standard ACDM model and the right panel for 
TeVeS. The dashed curves in the left panel show the line along which the integral is performed in (k, z) space for 
two values of I. Larger values of I probe smaller scales and therefore higher fc-values. At low z, the two theories 
make similar predictions (the larger amplitude for TeVeS at low redshift is model-dependent), whereas the differences 
become large for redshifts larger than one. 

The striking differences in between the two panels in Fig. 2 are due to the differences in the growth of structure 
in TeVeS and standard ACDM during the matter dominated era. In absence of cold dark matter, the growth factor 
in TeVeS must be enhanced with respect to ACDM in order to match the present amplitude of perturbations. This 
enhancement produces an evolution of the gravitational potential, and thus an ISW signal during matter domination. 
This does not happen in ACDM, where the ISW contribution comes only from the decay of the potential during dark 
energy domination. 

At redshifts of order one, Disw even becomes negative for TeVeS, leading to negative cross-power. For large 
z samples, this feature of TeVeS should have observable consequences. Before examining these consequences, we 
explain the features of TeVeS that lead to the behavior depicted in Fig. 2, in the following section. 
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FIG. 2: Left Panel, k 2 Djsw / Hq as a function of wavenumber k and redshift in standard ACDM. Here we show results from 
linear theory so the relevant growth functions do not depend on wavenumber. Dashed curves show the line along which the 
integral is performed for two values of I. Right panel, k 2 Djsw / Hq as a function of wavenumber k and redshift in TeVeS. Note 
that at high reshift, the TeVeS weighting function goes negative. 

III. THE GENERALIZED POISSON EQUATION IN TEVES 

Perturbations around the smooth cosmology in TeVeS were first studied in [10] with the full set of evolution 
equations derived in [11]. The equations which determine the scalar potentials of equation (1) are [11]: 



-2k 2 t> - 2e 4 ^ ( 3$ + fc 2 C + 3^* j = K B k 2 E + 8nGa 2 pS 



= ATiGa 2 e-^{p + P)Q 



(12) 
(13) 



where (j> is the zero order smooth part of the TeVeS scalar field; b = ae^; a and E characterize perturbations to the 
vector field; p is the background matter density; O is proportional to the velocity of the matter; 



C = (e 



-Acj> 



l)a; 



(14) 



and Kb is a dimensionless parameter governing the importance of the vector field. Here we have neglected the scalar 
field perturbations as these are small. 

In the following we will assume matter domination, including baryons and massive neutrinos (no cold dark matter). 
Using equation (13) we can write the term in parenthesis in equation (12) as: 



3$ + fc 2 C + 3-* j = 127rGa 2 e^(p + P)9 + k 2 (. 



(15) 



On sub-horizon scales we can neglect the O term. Then, using the definition of ( and | = - + <p, we obtain the 
generalized Poisson equation: 



$ = - 



3ff 2 {Q b S b + n„g„) K B E 
2ak 2 2 



(16) 



The first term on the right is the standard contribution from matter, but the last two are unique to TeVeS. Pertur- 
bations to the vector field affect the behavior of the gravitational potential, destroying the simple 1-1 relationship 
between potential and density implicit in general relativity. 

To get an equation for the other scalar potential, we use the relation [11]: 



* = $ + e 4 ^ 



C + 2 ( - + 20 ) C 



(17) 
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FIG. 3: Left Panel. |fc 2 "J>| in the TeVeS baryon model (see text). After integrating numerically the TeVeS perturbation equations 
we plot the output for |fc 2 $| and the different contributions to |fc 2 $| from the other perturbation variables, according to the 
Poisson equation. In the plot, "standard contrib." is the standard 8 term on the R.H.S. of the Poisson equation, "a contrib." 

represents the + ^ — e4 ^) a term and "E contrib." is the KbE/2 term (all these terms have been multiplied by k 2 ). 
"All terms" is the sum of all the R.H.S. terms. On the x-axis k is in units of h/Mpc but plotted lines correspond to units in 
which Ho ~ 1. Same for all the following plots. Right panel. Same curves, but in the TeVeS neutrino model. 



From equation (14) we get: 



The evolution equation for a is: 



hence: 



a = E + ty 



(e~^ - l) 



1 a. 



(18) 



(19) 



C = ~4e~ 4 *a + (e- A * - l) E + (e -4 * - 1 W + (V - I 



a 
a , 



We can now include these expressions into equation (17) to get, after some algebra: 



Using the Poisson equation (16) leads to: 

ie'^Hl (n b 5 b + n v s v ) 



- | e-^ - 1 
a 



) + 4> (e^ - 5 ) 



2ak 2 



: 1 llli | i 



E - 4(f>a 



(20) 



(21) 



(22) 



Again the first term on the right is familiar from general relativity but the vector field perturbations induce two new 
source terms. Summing $ and 'F, we obtain: 



$ + * 



3 Hi 



(fl b S b + nj u ) (1 + e-^) + 



Kf 



(1 + e" 4 *) + e 



4 ^ J_ - 4 * 



1 



E 



(- + 0)(1 - e-^) + . 
a 



(23) 
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FIG. 4: Left Panel. |fc 2 ^| in the TeVeS baryon model. After integrating numerically the TeVeS perturbation equations we 
plot the output for \k 2 ^f\ and the different contributions to \k 2 &\ from the other perturbation variables, according to the 
equation (22). Right panel. Same curves, but in the TeVeS neutrino model. 



The limit 4> <C 1 is imposed by nucleosynthesis bounds on fi^ [10]. 

Equation (23) can significantly differ from the Poisson equation in standard cosmology (equation (11)) if the vector 
field perturbations a and E are large. Two of us [12] showed that a and E get large when Kb —* 0. Thus we expect 
the Newtonian potentials $ and \T/ to deviate from standard Poisson expectations in this strong TeVeS regime. We 
also see that, unlike the standard term, the vector field terms in the TeVeS Poisson equation do not decrease as 1/fc 2 . 
This suggests that the corrections might actually dominate on small scales. 

We solved the full set of evolution equations in two different models. The first (called "baryon model" in the 
following), has only baryons with fib = 0.3 and Kb = 0.08; the second ("neutrino model") has massive neutrinos with 
Q u = 0.17, n b = 0.05, K B = 0.08. The left panel of Fig. 3 shows the contributions to |fc 2 $| given by the different 
terms in equation 16 for the baryon model. The right panel shows the same for the neutrino model. In both cases we 
can see that the standard term — 3e~ 4 ^_ff 2 (QmSm) /2ak 2 dominates on large scales but it becomes basically negligible 
on small scales, where the a contribution is the dominant one. An analogous result is obtained for |fc 2, F| (Fig. 4) but 
this time the dominant contribution on small scales is produced by the vector perturbation variable E instead of a. 

In Fig. 1 we showed the "Poisson ratio" — fc 2 (<f> + &)/ (f2&<5& + Sl v &v)- In a standard cosmological model this quantity 
is a constant but in TeVeS this is no longer true on small scales, where the vector field contributions in the Poisson 
equation dominate. In the baryon model the Poisson ratio shows large oscillations on small scales. The reason for this 
can be understood by looking at Fig. 3 and 4. When we consider the Poisson ratio on small scales we are basically 
dividing a and E (green and blue lines in the figures) by S ("standard contrib.", red line). We see that the acoustic 
oscillations in a and E are much more pronounced than and out of phase with those in 6. This produces the large 
oscillations that are finally observed. Acoustic oscillations are generated by the baryon component, which is much 
smaller in the neutrino model. This also explains why we do not see oscillations in the neutrino model Poisson ratio. 
Another feature of the neutrino model Poisson ratio is a decreasing trend on small scales (dashed-dotted blue line in 
Fig. 1). This is due to neutrinos entering the free streaming regime and thus producing a decay in S (red line in the 
right panels of Fig. 3 and 4) not matched by a corresponding decay in E and a. 

Finally, Fig. 5 (left panel) shows the effect of the vector field on the matter power spectrum. In this figure 
we considered our baryon model with Qb = 0.3 and we changed the value of Kb while keeping the other TeVeS 
parameters fixed as usual. We know that the effect of reducing Kb is to boost the growth of vector perturbations 
and to enhance the growth factor of matter perturbations as well (see [12]). We also see that reducing Kb shifts the 
peak of the power spectrum to smaller scales. 
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FIG. 5: Left panel: Matter power spectra for a TeVeS baryon only model with fixed = 0.3 and different values of Kb- Right 
panel: Galaxy-CMB cross correlation function at 9 = 6° (divided by the galaxy bias) from different surveys [30, 31]. The red 
(solid) curve shows the TeVeS neutrino model prediction, while the green (dashed) curve is the ACDM prediction. 

IV. RESULTS 

The authors of [10] showed that a flat TeVeS model with f2„ = 0.17, = 0.05, fl\ — 0.78 (i.e. the "neutrino 
model" we considered in the previous section) can produce a galaxy power spectrum in reasonable agreement with 
present observations. This means that a direct measurement of the power spectrum from galaxies is not the best 
discriminator of the TeVeS model and standard ACDM cosmology. It is then worthwhile to ask whether there are 
other large scale structure observables that can be used to distinguish between TeVeS and standard ACDM, even in 
the case of similar predicted power spectra. 

We know that in order to match a ACDM power spectrum without resorting to a CDM component, TeVeS needs 
a substantial contribution from the vector field perturbations E and a. As shown in the previous section, this leads 
to important modifications in the Poisson equation. These modifications are in turn responsible for the change in 
Djsw with respect to standard GR illustrated in Fig. 2 and can thus potentially produce an observable effect in the 
galaxy-CMB cross correlation given by equation (9) and (10). Note that since the TeVeS model contains a dark 
energy component, we also expect a late time ISW effect similar to that in the standard ACDM scenario when this 
component starts to dominate. 

For the following calculations, we use the ACDM linear power spectrum calculated using the transfer function from 
[32] for both TeVeS and ACDM, as it is accurate to the percent level, whereas the power spectrum in the TeVeS 
framework is not yet known to the desired degree of precision. Throughout, we assume the following parameters for 
the ACDM cosmology: h = 0.7, Q m = 0.27, Q A = 0.73, n b = 0.046, a 8 = 0.8, n s = 0.95. 

So far, the ISW effect has been detected in the cross-correlation of the CMB with galaxy surveys in the radio, 
infrared, optical, and X-ray bands (see [30] for a summary of recent results). These detections span a redshift range 
from 0.1 to 0.9. Recently, a cross-correlation with the SDSS quasar catalog has been performed [31] which has a 
median redshift of 1.5. Fig. 5 (right panel) shows the different observational results in terms of the angular cross- 
correlation function evaluated at 6 = 6°. Also shown are predictions from the TeVeS neutrino model and the ACDM 
cosmology defined above. In order to evaluate equation (10), we assumed a z-dependent galaxy selection function 
according to equation (12) in [30]. 

Clearly, the low- z (late-time) ISW effect predicted by TeVeS is quite similar to that of ACDM and consistent with 
the current observations. The effects of the TeVeS vector fields become noticeable only beyond z = 1, where there 
are no sufficiently strong constraints yet. 

In order to determine whether future surveys will be able to unambiguously distinguish the TeVeS scenario from 
ACDM, we adopt a fiducial survey ("Sample 1" in Ref. [29]) with parameters expected from the future Large Synoptic 
Survey Telescope (LSST, [33]): we assume a galaxy sample with a limiting magnitude of 27 in the I band and a sky 
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FIG. 6: Left panel: The cross-correlation of galaxies and the CMB for a LSST-like galaxy survey and four redshift bins, in 
TeVeS (thick lines) and GR (ACDM, thin lines). Right panel: The significance of the TeVeS — GR deviation in the galaxy-CMB 
cross power for the LSST-like survey; i.e. difference between TeVeS and GR divided by the expected error in each angular bin. 

coverage of f s — 0.5. Dividing the survey into several redshift bins, the observed number of galaxies per arcmin 2 as 
a function of redshift for each bin, W(z), is then given by the observed z-dependent luminosity function of galaxies 
[34, 35] convolved with a smoothed top hat redshift window (Eq. (16) in [29]). This roughly takes into account the 
scatter expected in the photometric redshifts. The mean redshifts of galaxies in the four bins considered range from 
0.49 up to 3.54. In order to convert from galaxy overdensity to matter overdensity, we have to assume a mean galaxy 
bias in each redshift bin. This was calculated in [29] using the halo occupation distribution and the halo mass function 
from simulations [36, 37]. In actual surveys, the bias is determined from the data itself. The characteristics of the 
redshift bins are summarized in Tab. I. With these parameters, we can evaluate equation (9). Note that we neglect the 
magnification bias [29] as well as redshift errors here. These effects are not expected to affect our main conclusions. 

Fig. 6 (left panel) shows the cross-power spectra for TeVeS (thick lines) and ACDM (thin lines) for the different 
redshift bins. As expected, the TeVeS cross-power is similar to ACDM in the lowest redshift bin and goes negative at 
high z (the "baryon model" considered above shows the same qualitative behavior). To estimate the detectability of 
this effect, we bin the cross-power spectra in /, taking the errors on a given multipolc moment as: 

(AC 3T ) 2 (Z) = /a(2 * +1) (c 2 gT (l) + C TT (l)(C gg (l) + -L)) , (24) 

where C g r is the galaxy-temperature cross-correlation; Ctt and C gg are the CMB and galaxy auto-power spectra 
respectively; and n g is the number density of galaxies per sr (so that 4irf s n g is the total number of galaxies in the 
corresponding redshift bin of the survey; see Tab. I). Throughout we set f s = 0.5, and errors on different C 9 t{1) are 
assumed to be uncorrelated. Fig. 6 (left panel) then shows the difference between the adopted ACDM model and 
TeVeS divided by the expected error in each angular bin. The TeVeS predictions are clearly distinguishable from 
the ACDM model, reaching significances of over 10cx in the high redshift bins. In addition, the negative sign of the 
correlation at high z is an unambiguous signature of this modified gravity theory. 





Bin 1 


Bin 2 


Bin 3 


Bin 4 
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0.49 


1.93 


2.74 


3.54 
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1.08 


2.02 


2.90 


3.89 


n g 


34.6 


10.1 


3.89 


1.68 



TABLE I: Parameters of the four redshift bins for the assumed LSST-like survey [29]. (z) is the average redshift of galaxies in 
the bin, b is the mean galaxy bias, and n g is the observed number of galaxies per arcmin 2 in each redshift bin. 
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V. SUMMARY AND CONCLUSIONS 



TeVeS is a modified gravity theory originally introduced by Bekenstein [9] with the purpose of providing a covariant 
relativistic framework for Milgrom's paradigm of MOdificd Newtonian Dynamics (MOND). Further investigations [10] 
showed that a TeVeS model with massive neutrinos and no cold dark matter component can pass some fundamental 
cosmological tests. In particular, it was shown that such a model can reproduce observations of the CMB and galaxy 
power spectra. 

In light of these results, in this paper we raised the following question: are there cosmological observables that are 
able to clearly distinguish between TeVeS and standard ACDM ? 

In order to reproduce a ACDM power spectrum without using cold dark matter, TeVeS has to incorporate a 
mechanism to enhance the standard growth rate of perturbations [12]. This in turn produces a characteristic ISW 
signature detectable in principle by looking at the CMB-galaxy cross-correlation. In section II we computed the ISW 
effect in TeVeS, showing that the characteristic quantity Djsw(k, z) presents clear differences between TeVeS and 
ACDM at intermediate redshifts z > 1. By deriving the generalized Poisson equation in the TeVeS framework (Section 
III) we showed that the same vector field perturbations that are responsible for the growth of large scale structure in 
TeVeS are also the cause for the different predictions for the ISW effect. 

The TeVeS model studied here is consistent with current observations of the ISW effect (Fig. 5, right panel) which 
reach up to z w 1.5. In order to determine whether the TeVeS predictions are constrainablc in future surveys, we 
computed the expected CMB-galaxy cross-correlation for a fiducial survey with parameters expected from the future 
LSST (Section IV). Our results show that at high redshifts [z > 2), TeVeS and ACDM will be clearly distinguishable, 
with significances of over 10a. Moreover we found that the sign of the CMB-galaxy cross-correlation at high redshifts 
differs between TeVeS and ACDM, thus providing a smoking-gun signature for this kind of modified gravity. We 
emphasize that this CMB-galaxy anti- correlation at high redshifts is a robust prediction, since it is due to the 
additional degrees of freedom that are necessary for large scale structure to form in a TeVeS cosmology. 

We would finally like to stress that our work was conducted in the framework of TeVeS because this provides a 
well-defined test model for which the evolution of linear perturbations has already been studied. However, as pointed 
out by Bertschinger [20], the characteristic TeVeS features that give rise to the studied effect are expected to be 
common to a large class of modified gravity models. Thus we expect the CMB-galaxy cross-correlation to be a general 
and powerful tool to discriminate among many alternative modified gravity theories. 
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